Global spectral representations of black hole spacetimes in the complex plane 
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Binary black hole coalescence produces a finite burst of gravitational radiation which propagates 
towards quiescent infinity. These spacetimes are analytic about infinity and contain a dimension- 
less coupling constant M/s, where M denotes the total mass-energy and s an imaginary distance. 
This introduces globally convergent Fourier series on a complex radial coordinate, allowing spectral 
representation of black hole spacetimes in all three dimensions. We illustrate this representation 
theory on a Fourier- Legendre expansion of Boyer-Lindquist initial data and a scalar wave equation 
with signal recovery by Cauchy's integral formula. 
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I. INTRODUCTION 

Rapid progress in sensitivity in the broadband 
gravitational-wave observatories is creating new oppor- 
tunities for searches for the bursts of gravitational-wave 
emissions produced by coalescing black holes. Detection 
strategies will benefit greatly from a priori understand- 
ing of their gravitational-wave emissions, which includes 
matched filtering techniques against a catalogue of sig- 
nals precomputed by numerical relativity. This poses an 
interesting challenge of designing highly efficient and sta- 
ble computational algorithms. 

Numerical relativity is particularly challenging in cal- 
culating bursts and preceding chirps over many wave- 
periods in the presence of singularities associated with 
black hole spacetimes. Spectral methods provide an at- 
tractive approach. They are optimal in efficiency and ac- 
curacy, both in amplitude and phase, provided that the 
metric is smooth everywhere. Here, we focus on bursts 
of radiation produced by black hole coalescence which 
propagate towards quiescent infinity. By causality, these 
spacetimes preserve quiescence at infinity for all finite 
time. 

The asymptotic structure of black hole spacetimes 
which are asymptotically flat and quiescent at infinity 
shares the same asymptotic properties as the Green's 
function of Minkowski spacetime. Consequently, these 
spacetimes carry a dimensionless coupling constant M/s, 
where M denotes the total mass-energy of the space. The 
Green's function of the Laplacean, 
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is a function of the distance D between an observer at y l 
and a source at p % . We can expand G in Legendre poly- 
nomials Pi using a spherical coordinate system (r, 9, </>), 



G = 



yfr 2 + p 2 — 2pr cos 9 



which shows that G is analytic in each coordinate away 
from the singularity at y l — p l with domain of conver- 
gence r > p. We introduce complex distances through a 



complex radial coordinate 



z = x + is (s > p). 



The analytic continuation 
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G = p- 1 E(^) P 2i (cos( 
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is globally convergent on — oo < x < oo outside \z\ = p. 
(This corresponding branch cut for the square root is 
(—oo,0].) To be precise: for each 9 Q defines a Taylor 
series in 1/z and, for each z, it defines a convergent ex- 
pansion in Legendre functions. Thus, p/s is a running 
coupling constant and it must be less than one for con- 
vergence. We shall refer to M/s as the fixed coupling 
constant, noting that p/s is decreasing in binary coales- 
cences. 

Analyticity of the metric in lm(z)> s defines an op- 
timal distribution of points on z = x + is for numerical 
implementation by the conformal transformation 



to 



M 

z 



(5) 



Here, w describes a circle of radius (2s) _1 about wq — 
(2is)~ 1 , whose interior corresponds to Im(z)> s. There 
exists a Taylor series in w about wq which is convergent 
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up to including the boundary C : w 
The corresponding Fourier series in (f> can be efficiently 
computed using the fast Fourier transform on a uniform 
distribution of points (j) k = 2- m+1 Trk {k = 1, 2, • • • , 2 m ) 
on C. Their images on z = x + is produce the non- 
uniform distribution 



z k = Ms [tan(0 fc /2) + i] , 



(6) 



whose large-scale cut-off is set by the proximity of the 
{> fc } to ±7T. 

We apply these observations to Boyer-Lindquist initial 
data for two black holes 0, , described by the potential 
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\J z 2 + p 2 — 2pz cos 6 \J z 2 + q 2 — 2qz cos 9 
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for a pair of black holes with masses Mi t 2 at positions 
p and q along the polar axis of a spherical coordinate 
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FIG. 1: The Boyer-Lindquist scalar field &(z,8) possesses a Fourier-Legendre expansion in w = l/z,z = x + is, (top two 
windows) and cos 8 (bottom two windows), here shown for a symmetric configuration consisting of two black holes of mass 
Mi,2 = 1 positioned at the north and south pole p = ±1 of a spherical coordinate system (r,8,(j>). The Fourier coefficients 
show exponential decay by analyticity in Im(z) > s > s*(0) — psind shown with cos# = 0.25 and s* = 0.9682. This gives rise 
to spectral accuracy in the first and second coordinate derivatives (top right). The norm of the Legendre coefficients shows 
exponential decay when s is larger than p (left). This gives rise to spectral accuracy in the first and second derivatives with 
respect to 8 (bottom right). 



system (r,8,<fi). It will be appreciated that the singu- 
larities z — pe ±l6 are essential in making $ non-trivial: 
the alternative of an entire function which is analytic at 
infinity would be constant by Liouville's theorem. For 
complex z = x + is with s larger than both p and q, 
<fr(z,9) is finite and pointwise analytic everywhere on 
— oo + is < z < oo + is. Following Q and ©, we capture 
analyticity by the Fourier-Legendre transform on the im- 
age w = l/z. This gives rise to an efficient and accurate 
functional representation in (z, 9) with spectral accuracy 
in derivatives by spectral differentiation. Fig. 1 shows 
the computational results on the alternative representa- 
tions 

$(z,0) = Z% =1 a n (cosd)w n - £f =o Mz)P ; (cos0), (8) 
where the series are finite in our numerical implementa- 



tions. 

The application to black holes spacetimes in the com- 
plex plane can be considered in the line-element 



-N 2 dt 2 +h i] dx i dx j , N(z) 



1 - M/2z 
1 + M/2z ' 



(9) 



for a three-metric h%j. Working in the complex plane, 
we are at liberty to choose a lapse function N and 
vanishing shift functions in accord with the asymptotic 
Schwarzschild structure of spacetime at large distances. 
This applies without loss of generality to any binary black 
hole spacetime with total mass-energy M as measured at 
infinity, defined by the residue 
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FIG. 2: Normalized errors in the first (circles) and second (squares) derivatives of the metric and the Ricci scalar tensor (stars) 
based on a Fourier-Legendre representation of Boyer-Lindquist initial data of Fig. 1 with M/s — 1/2. The Fourier-Legendre 
results are computed with N = 32 Fourier coefficients on w = 1/z varying degrees L = 4, • • • ,32 of Legendre polynomials. 
The results are compared with a mized finite differencing- Legendre method with varying degree N = 2 m , m = 5, • ■ ■ 11 of 
the number of in the z— coordinate and fixed expansion in L — 16 Legendre polynomials. This comparison highlights the the 
accuracy and efficiency of the Fourier-spectral method, even at low resolutions. The results are practically unchanged for other 
configurations of Boyer-Lindquist data, and the errors decrease for larger s. 



For binary black hole coalescence, M is the sum Mi + A/2 
of the masses of the two black hole, plus the energy in 
gravitational radiation and minus the binding energy in 
the system. Thus, M is a constant. In contrast, M is 
a non-increasing function of time on any finite computa- 
tional domain with outgoing radiation boundary condi- 
tions. For large z, the line-element (0 satisfies 



ds 2 = -N 2 dt 2 



1 + M.) (dz 2 + dY?) 
2z I 



(11) 



up to order (M/z) 2 , where dT? = z 2 d9 2 
notes the line-element of the coordinate sphere of radius 
z. Accordingly, the three- metric hij is decomposed into 
diagonal and off-diagonal elements given by the matrix 
factorization 



A\ B x B 2 
h tJ = 7 ( B, A\ B 3 1 7) 
B 2 B 3 Af 

where 7 is the diagonal matrix 

711 = 1,722 = 2,733 = zsin( 



(12) 



(13) 



In accord with the spectral representation of @ and the 
potential Q, the coefficients (Aj,Sj) are expanded in 
spherical harmonics on z = x + is, 



Ai = M£| m |<; af(t)w l+1 Y ml ( 



B, 



M£ M < ; bf(t)w l+1 Y n , 
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Together with the Fourier expansion onui = M/z, 1141 
1151) defines a spectral representation in all three coordi- 
nates. Fig. @ shows the results in case of the axisym- 
metric Boyer-Lindquist initial data, including the compu- 
tational error behavior in the scalar Ricci tensor. Here, 
we use the general expansions (JSJ) to probe independently 
the dependencies on N and L (to go beyond N = L + 1 
in any truncated series of H14I15H L 

Time-evolution of hij on z = x + is gives rise to burst 
of radiation which propagates towards quiescent infinity, 
where the observer's signal is defined on non-negative 
radii z > 0. This signal can be calculated by projection 
according to Cauchy's integral formula, 



u(r, 9, 0, t) = Re 



oo+is 



u(z, 9, 4>, t) 



dz 



— oc-Hs 



(16) 



where u denotes a relevant metric component or tensor 
quantity. Time-evolution on z — x + is is stable when 
constant translation s, since this preserves reality of the 
dispersion relation. This can be illustrated by the linear 
wave-equation for a scalar field u(t, r) in spherical coor- 
dinates (t, r, 9, <fi) with time-symmetric initial data, given 
by 



K 



u tt {t,r) =r 2 [r 2 u r (t,r)] r 
u(0,r) = u (r), u t (0,r) = 0. 



(17) 



By analytic continuation, the equivalent initial value 
problem K' with complex radial coordinate on — 00 +is < 
z < 00 + is is 



K 



u tt (t,z) ^w A u ww (w = l/z) 
u(0,z) =uq(z), Ut(0,z)=0. 



(18) 
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FIG. 3: Signal recovery of a wave propagating to quiescent infinity is shown for a quantity u(t, z) on z = x + is, here satisfying 
the scalar wave equation u« = r~ 2 (r 2 u r ) r with time-symmetric initial data at t = u(0,x) = e~ 2x cos(fCte) and Ut(0, x) — 0. 
The simulations use the Fourier representation on the complex radius under w = M/z and the wave-equation utt = w u ww 
using N = 1024 points. Plotted is the quantity ru on r > for to — and t\ — 4. The numerical results (dots) accurately track 
the exact results (continuous line, by the method of characteristics). The errors are essentially independent of s, i.e., 4.5e-4 for 
s = 0.5 and 4.7e-4 for s — 1.0, as implied by analyticity in s > 0. 




(In general, the analytic continuation of the initial data is 
the extension of data on (r, 9, cf>) and (r, tt~6, </>+7r)). We 
can implement K' numerically using finite-differencing 
and leapfrog time-stepping. Fig. © shows a represen- 
tative numerical result, which is verified to satisfy in- 
dependence of s. The scalar wave-equation with time- 
symmetric initial data shows the appearance of both 
left and right movers, the first propagating towards z = 
— oo+is and the second propagating towards z = oo + is. 
Only the right movers are projected onto the observer's 
non-negative radii z > 0. The left movers are ultimately 
dissipated "unseen" in their propagation to — oo + is, 
while maintaining a contribution to the total mass-energy 
M in CT - 

To conclude, the complex plane offers a unique op- 
portunity for circumventing spacetime singularities and 



enabling spectral representations in all three dimensions, 
consisting of spherical harmonics in the angular coordi- 
nates and Fourier series in the complex radial coordinate. 
This approach creates the most efficient representation 
of the metric, promising a reduction in computational ef- 
fort by orders of magnitude compared to standard finite 
differencing approaches. It will be of interest to take ad- 
vantage of this approach in time-dependent calculations 
described by the fully nonlinear equations for general rel- 
ativity in any of its hyperbolic forms (e.g., 0)- 
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